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We have analyzed the interplay between noise and periodic modulations in a mean field model of 
a neural excitable medium. To this purpose, we have considered two types of modulations; namely, 
variations of the resistance and oscillations of the threshold. In both cases, stochastic resonance is 
present, irrespective of if the system is monostable or bistable. 
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I. INTRODUCTION 

When a nonlinear system is driven by a periodic force in a noisy environment, its response may be enhanced by the 
presence of noise. This constructive role played by noise can be characterized by the appearance of a maximum in the 
so called signal-to-noise ratio (SNR) at a nonzero noise level. In essence, the SNR is a quantity that reflects the quality 
of the output signal, in such a way that for large values of this quantity the output signal looks more ordered. This 
phenomenon, named stochastic resonance (SR) has been found in many situations pertaining to different scientific 
areas |^2(J. 

In regards to neural systems, many examples of SR has been found theoretically in single neurons J2l],^2| and neural 
networks |2^-|2^], and experimentally in single neurons [ pGp^l. In contrast, an aspect which has only been considered 
recently is the appearance of SR in neural excitable media |28| . In a set of experiments with mammalian brain slices 
corresponding to the hippocampal center ]29[ ] it was demonstrated that an electric field can either suppress or enhance 
coherent activity in real networks. It was later shown, by using a time- varying electric field, that as the magnitude of 
the stochastic component of the field was increased, SR was observed in the response of the neural network to a weak 
periodic signal [^8| . These and other recent results ||(| clearly shows that neural noise could play a relevant role in 
the information processing of the brain. 

In this context, both key ingredients for SR, i.e. a well-defined coherent and time-periodic modulation and intrinsic 
noise, are present at several scales in neural tissues. We can consider a coarse-grained characterization of brain 
dynamics by means of the analysis of the electroencephalogram, which is an averaged measure of the spatiotemporal 
activity of millions of neurons. These neurons, in turn, are part of a network receiving inputs from various parts of 
the nervous system called nuclei. Among these, the thalamus plays a very important role in controlling the behavioral 
states of the brain |3l| . In fact, the thalamus is known to display autonomous oscillations p]J| , i.e. it works as some 
kind of pacemaker to the brain cortex. The thalamus is massively connected with the cortex and produces autonomous 
periodic oscillations, even when disconnected from the cortex. Hence, the cortex can be, to some extent, viewed as a 
highly connected network which is periodically stimulated by the thalamic pacemaker. Several experimental results 
give support to the idea that the behavioral states of the brain (alpha rhythm, slow delta waves, sleep, and REM) 
are somehow related to an oscillatory input into the cortical tissue |32| . These facts indicate that a realistic model of 
cortical dynamics should consider the effects of the thalamo-cortical pacemaker, which can be simulated in different 
ways, but it can always be considered as a periodic external signal. 

Another aspect to be emphasized is the fact that single neurons of the cortex exhibit some degree of variability, 
i.e. the response to a stimulus depends on the particular trial, which is observed in single-neuron experiments and 
also in measurements of single-neurons inside the brain. Such a variability comes from both complex deterministic 
dynamics and the noise implicit in the random nature of the incoming signals. Additionally, single neurons manifest 
intrinsic noise. In this sense, there have been numerous experimental studies about the stochastic activity of nerve 
cells. Noise has been observed in nerve-cell preparations and involves both synaptic noise, which is due to randomly 
occurring synaptic potentials p3), and membrane noise. In the last case, small fluctuations in the electric potential 
across the nerve-cell membrane are observed, even when apparently steady conditions prevail. These fluctuations are 
linked with conductance changes induced by random closing and opening of ion channels. 

The aim of this paper is to study a simple model which is able to capture the main traits of the actual neural media, 
as concerns the interplay between noise and periodic modulation in neural dynamics. To this purpose, we present a 
detailed analysis of a standard mean field model of neural excitable medium which was introduced in a preliminary 
form in Ref. @. 
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The paper has been organized as follows: In section |l] we present the mean field model. In sections III and [V, we 
analyze the appearance of SR in two different situations concerning oscillations of the resistance of the neural tissue 
and variation of the threshold. Finally, in section ^ we summarize our main results and discuss possible implications 
in neural systems. 

II. MEAN FIELD MODEL 

Let us consider a slab of neural tissue comprising a very large number of closely packed and coupled nerve cells, where 
connections are only excitatory. Different parts of the neocortex can fit this description with more or less success. In 
general, most of the real networks formed in the brain cortex are constituted by both excitatory and inhibitory neurons. 
The frequency and relevance of each type of cell, however, vary depending on the place considered. In this sense, 
the overall activity of some brain regions, like the thalamus, is largely dominated by the excitatory component |p4) , 
usually identified, at the microscopic level, with the piramidal cells. This excitatory behavior is particularly relevant 
in the hippocampus, when stimulation of a single excitatory neuron can generate a burst of synchronous activity |p5ft . 
These results makes the consideration of the behavior of cortical nets in terms of purely excitable dynamics reasonable, 
as a first approximation. In addition, the lack of a microscopic structure is not an important problem when neurons 
and connections are not explicitly considered as discrete entities, as done at our level of description. In this regard, 
previous studies reveal that the mean-field approach is able to reproduce the detailed macroscopic description [p6[ . 

The situation we consider involves one of the simplest models of aggregates of nerve cells. At a given spatial point r, 
the quantity of interest is the mean local potential V{r) which is the result of a local integration of incoming signals, 
i. e. 

V(f,t)= J v{f-f',t)p(r')df' , (1) 

over a given neighborhood T. Here p(f) is the local packing density of neurons and v(f, t) the transmenbrane potential. 
The time evolution of V(r, t) can be obtained from the following integrodifferential equation 

c BV^ = _^ +mHvm) + im (2) 

Here C is the capacitance, R the resistance, I(r,t) a stimulus and (3{r) the mean number of synaptic connections. 
The functional form of $ is defined through the sigmoidal relation: 

$(V) = [l + e-^-^]" 1 , (3) 

where $o an d v are given constants and is a threshold. Many possible functional forms for j3(r) can be considered, 
as for example an exponential decay j3(r) = 6exp(— |r|/7) which leads to the equation: 

dV(f,t) _ V(f,t) 
dt R 



l + e 



-v{V(r' ,t)-6) 



dr' 



+ I(r,t) . (4) 

We will consider the case of all-to-all connectivity, i. e. 7 L, where L is the characteristic length scale. Under 
this assumption, the model of neural excitable medium constitutes a mean field approximation and is formulated by 
the equation |37| 

C^- = -R-'V + e{\ + e-'^-^)- 1 + I , (5) 

accounting for the dynamics of the spatial average V of the transmembrane potential. Here e is a constant that arises 
from that approximation. This model is referred to as the Cowan-Ermentrout model (CEM). Its dynamics can be 
described through a potential function U(V), 

dV dU I 

iir = -w + c- (6) 
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In our case, 



U(V) = — vR- 1 V 2 -2e\n [l + e v{ y- e A . (7) 

A remarkable fact is that for some values of the parameters, CEM exhibits bistability. Therefore, small changes in 
the parameter values can lead to sudden shifts from one stable branch to the other |p7| . 

To render our analysis of the mean field model complete, we need to specify the nature of the noise. It is worth 
pointing out that its origin, its characterization, and its effects in actual neural media inside the brain is far from being 
clear. Here, in order to account for the noise effects, we will consider a simplified situation that can be described by 
a fluctuating current applied to the net. Moreover, we assume that this current may be approximated by a Gaussian 
white noise [(/(*)) = and (I{t)I(t + r)) = 2a6(r)]. 

Another remarkable fact is that periodic modulations may affect the system in different ways, for instance, by 
periodically changing the value of one of the parameters. Thus, small changes in the permeability of a suitable ion 
gives rise to variations of the membrane resistance R. Electric fields may cause shifts in the effective threshold 9 for 
an action potential initiation. The parameter 6 could also change when stimulus are sent to the medium from other 
regions, e.g. when two networks interact. Here we will explicitly consider these two cases, i.e. oscillations of the 
resistance and variations of the threshold, although other possibilities could also be analyzed. 

III. OSCILLATIONS OF THE RESISTANCE. 

As a source of periodic modulation, we will first focus on the oscillations of the membrane resistance. In order to 
proceed with our analysis, we will consider the Fokker-Planck equation giving the probability density of the spatial 
averaged transmenbrane potential. For the sake of simplicity, we introduce the new set of variables: x = V, Xq = 0, 
k(t) = R(t)~ 1 C~ 1 , D = <t/C, and e = e/C. In these variables, the Fokker-Planck equation reads 

dP d d 2 P 

-=-- l{Fo + cFl)P] + D —, (8) 

where Fq — —k(t)x, with k(t) = k[1 + acos(woi)], is the linear force and F\ = [1 + e~ v ^ x ~ Xo >]~ 1 accounts for the 
nonlinear contribution, k, a, and ujq are constant parameters. 

The probability density can be expanded in powers of the strength of the nonlinear term; namely, P = Pq + eP\ + 
£ 2 Pi + . . .. Here Pq corresponds to the linearized equation whereas the remaining terms account for corrections 
due to the nonlinearity. Substituting this expansion in Eq. (|8|) we obtain the evolution equations for the different 
contributions: 

dP °~ d (FP) + D d2Po fQl 

~dT--dx- {FoPo) + D ^ ' (9) 

!> ' } ' 3 (F P 1 +F 1 P )+D^ , (10) 



at ox ox 1 
dP 2 _ d_ 

at ox ox 1 



-^ = -in + F1P1) + D—i . (11) 



To proceed further, we will define the quantities 

(x t ) i = I xP t (x,t)dx , (12) 



with i = l,2, 3, and 

/oo />oo 
xdx / yP i (y,t + T\x,t)P 3 (x,t)dy , (13) 
-00 J —oc 

with i 7 j = 1,2,3. From Eqs. (|9|)-([Tl|) one can easily see that (xt) = 0, whereas (xt) 1 ^ and (xt) 2 7^ 0. Moreover, 
the correlation function can be expanded in the form 
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xdx I yP(y,t + T\x,t)P(x,t)dy 
-OO <J — OO 

+ e (x t+T xt) 10 + e (x t+T x t ) l 

+ £ 2 { x t+rXt)i i 

+ e 2 (x t+T x t ) 2 Q + e 2 {x t+ rXt) Q 2 ■ ( 14 ) 
To analyze the interplay between noise and input signal, we will consider the SNR, defined as usual by 

where the output noise is in first approximation given by 



(16) 



N(oj) = — / / (x t+T x t ) cos(uiT)dTdt , 

w Jo J~oo 

and the output signal S(luo) comes from 

S(uj ) [5(uq -cj) + 6(ujo + uj)} = 
= — / / e 2 {x t+T ) 1 {xt) 1 cos(ujT)dTdt . (17) 

^0 Jo J-oo 

Since for large r the quantities Xt+ T and Xt become uncorrelated, in the previous equation we have replaced {xt+rXtjj^ j 
by {x t+T ) i (xt)j due to the fact that to compute the signal it is sufficient to know the behavior of the correlation function 
for larger times. Therefore, to obtain the SNR it is sufficient to consider only Po(y, t + t\x, t), Pq{x, t) and Pi (a;, t). 

With the purpose of obtaining the equation for the probability density, we will assume that the resistance varies 
slowly and that the amplitude of the oscillations is small. By using these approximations and by taking into account 
the fact that when the contribution proportional to a can be neglected the linear system constitutes an Ornstein- 
Uhlenbeck process, the noise term can be readily computed from 

(x t+T x t ) Q0 = -e~ KT . (18) 



Therefore the spectral density of the output noise is 



VU)= 2^' (19) 



provided that uj <C k. 

Moreover, we will also assume that the sigmoidal function [Eq. (|^)], giving the mean firing rate, may be approx- 
imated by a step function. This fact occurs when the gain of the neuron v is sufficiently large. In such a case, the 
potential function is given by 

U(V) = ^ [R-'V 2 - 2e(x - x Q )Q(x - x a )] . (20) 

Under these circumstances, 

p _ ]_ e -k(t)x 2 /2D e E{x~x )e(x-x )/D ^21) 

where Z is the normalization factor. Up to order e Eq. ( pl| ) yields 

P = P +sP 1 , (22) 

where 



l] < = \I^D e ~ mX2/2D ' (23) 
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and 



ft = iPo 



(a; — xo)0(x — xq) — / (x — Xo)Podx 



Notice that by using the fact that a is small Po can be approximated by 

Po = 



-kx 2 /2D 



2ttD 



1 



1 + |(1 - x 2 K,/D)acos(uj t) 



In order to obtain the output signal, we will take into account the expression 

{x t ) 1 = B + Aa cos(woi) j 



(24) 



(25) 



(26) 



which holds when a and loq are sufficiently small. Here B and the A do not depend on time. By using this expression 
in Eq. (^) we obtain the signal as a function of the susceptibility A, 



This quantity can be computed from Eqs. fl25| ) and (|2J), and one obtains 

A 



1 



«_ -kx*/2D , x _ \n_ x 2 K l D \ dx 

2D V 2ttD v on / ) 



The SNR then reads 



SNR = n(aeA) 2 — 



(27) 



(28) 



(29) 



The previous integral cannot be performed explicitly, however, its behavior for high and low noise levels can be 
obtained easily. 

The high noise level case can be performed by replacing the lower limit of the integral xq by 0, provided that 
Xq <C D/k. We then obtain 



a = — x p^ r z- kx2/2d (* 2 - d*=-- 

2D V 2ttD J V D J n 



Therefore 



and 



2 K 



SNR = 7r(ae) s 



D 



In the same way, for low noise level we can perform an asymptotic expansion by using the formula 



i n e~ au du= -x n+1 - 



ax' 



In this case the susceptibility is given by 



A = —x Q 



1 



2-kkD 



71-1 

2ax 2 



-kxI/2D 



Therefore, the signal and the SNR are 



S(uj ) 



(x ae) 2 e~ KX ^ D 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



and 
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WR= l -(x^)\e-^' D . (36) 

At low noise level the SNR increases whereas it decreases for large values of the noise, therefore the SNR has a 
maximum which indicates the presence of SR. 

In order to analyze the case in which the oscillations and the nonlinear term are small, but not infinitesimal, we 
have numerically integrated the corresponding equations following a standard second order Runge-Kutta method for 
stochastic differential equations ]3^ , |39t . The Langevin equation we have integrated is the one that corresponds to the 
Fokker-Plank equation [Eq. (||)] and is given by 

J = -k(t)x + e[l + e -"(— «)]-! + £(t) , (37) 

where £(t) is Gaussian white noise with zero mean, and correlation function = 2D6(t — t'). Here, small 

means that the effects of the nonlinear term are not large enough in order for the system to become bistable. The 
potential function giving the dynamics of the system is depicted in Fig. |l|a, for different values of the resistance. In 
Fig. ||a we have shown the behavior of the SNR as a function of the noise level D for two frequencies. The values 
of the remaining parameters are the same as the corresponding ones to the potential function of Fig. |l|a. This figure 
clearly exhibits a maximum in the SNR. 

When the nonlinear term is large enough, the system may display bistability. One state corresponds to all neurons 
at rest and the other to active neurons. 

In Fig. |3|a we have displayed the potential function associated to the Eq. (|3?]), when the resistance varies for 
values of the parameters corresponding to the bistable situation. Note that when the resistance depends on time, the 
position of the minimum corresponding to the active state also changes periodically in time. 

In order to analyze this situation, we have numerically integrated the corresponding equations as in the previous 
situation where the potential function is monostable. Fig. |4|a displays the SNR which exhibits a maximum at a 
nonzero noise level. 

In view to illustrate how the system behaves, in Fig. ||a we have shown three time series for different noise levels. 
This figure clearly manifests the presence of an optimum noise level, at which the response of the system is enhanced, 
and the displacement of the minima corresponding to active neurons. 



IV. VARIATION OF THE THRESHOLD 



In the previous analysis we have studied the case in which the resistance of the neuron undergoes oscillations. 
Another possibility of temporal variation are oscillations in the parameter 9. Explicitly, the dynamics corresponding 
to this situation is described again by Eq. (g), but now F = —nx and Fx = [1 + e - t/ ( 2; - a: o- QCOS ("ot))]-i^ gy us j n g 
the same assumptions about the parameters a, e, ujq, and the gain v introduced previously, we will now proceed in a 
similar fashion. 

It is easy to see that in this case Pq does not depend on time: 



Po = ^2^D e D ' (38) 



and the correction to the probability density due to the nonlinear term is given by 

Pi = 7Tp P o{[ x -xq- acos(uj Q t)]Q(x - x - acos(ui t)) 



[x — xq — a cos(u>ot)]Po(x)dx} . (39) 



The averaged value of x is then 



f°° 1 

(x) 1 = / x— P (x)[x - x - a cos(u> t)]dx . (40) 

J xo+a cos(cjot) 

Note that for symmetry reasons, the integral in Eq. (B3) gives a null contribution to (x) 1 . By expanding in the 
parameter a around its zero value, we then obtain Eq. (Eq). In this case the susceptibility reads 
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r°° i 1 

A = -\ x w> p ^ dx = —^=^ e ~ KxV2D ■ ( 41 ) 

The noise term is the same as in the previous case, then we obtain the following expressions for the signal, noise 
and SNR 

S{ij )^^(aeA) 2 = —^—(ae) 2 e- KX ^ D , (42) 



2 x ' WitkD 
D 
2^2 



iVM = — 2 , (43) 



SNR = n(aeA) 2 ^ = ^^e^/D . 



(44) 



This last expression clearly shows that the SNR has a maximum at a non-zero noise level then making the presence 
of SR manifest. 

As we did in the previous situation, to study the case in which the oscillations and the nonlinear term are not 
infinitesimal, we have numerically solved the corresponding Langevin equation by using the procedure outlined above. 
In this case it reads 

^ = -Kx + e[l + e- l,{x - Xa - acos{uot)) ]- 1 +£(*) , (45) 

where is the same noise as the one defined previously. The explicit situation we have considered is given through 
the potential function displayed in Fig. [l]b. In Fig. |^b we have shown the behavior of the SNR as a function of the 
noise level D for two frequencies. The values of the remaining parameters are the same as those corresponding to the 
potential function of Fig. |l|b. In this case the SNR also exhibits a maximum for the two frequencies. 

Bistability may also be present in this case. In Fig. |^b we have represented the potential function when periodic 
modulations act through the threshold, for values of the parameters corresponding to the bistable situation. In 
contrast with the case of oscillations of the resistance, in which the minimum corresponding to the active state varies 
its position, when the threshold oscillates, the two minima always remain at the same transmenbrane potential. 

This situation can also be analyzed through numerical integration. Fig. |]b displays the SNR for the periodic 
modulation we are considering. This quantity exhibits a maximum for the two frequencies. Finally in Fig |^b we 
have also shown three time series for different noise levels. It is worth emphasizing that, for noise levels close to the 
optimum value, the time series look as those corresponding to the usual bistable quartic potential ||[|. 



V. DISCUSSION 



In this paper we have analyzed how noise affects the behavior of a neural medium when it is periodically modulated. 
We have found that the occurrence of noise may play a constructive role since an optimized amount of it may 
contribute to enhancing the response of the system. Under some circumstances, the presence of noise is responsible 
for the appearance of oscillations which otherwise would not be manifested. In this regard, the analysis of macroscopic 
neural dynamics obtained from electroencephalograms has been a matter of debate over the last decade fi(| . It is 
accepted that the activity of the brain cortex shows low-dimensional traits, though the exact nature of the phenomenon 
itself is far from being clear. Here we have shown that noise, sometimes not considered, could give rise to a coherent 
behavior of the system, then playing an important role in neural dynamics. 

The mean field model we have proposed, although constituting an oversimplified picture of a thalamo-cortical 
network, might be a first step in our understanding of how noise and nonlinearities can generate interesting macroscopic 
outcomes. Further developments must include spatial effects as well as the consideration of activatory and inhibitory 
populations of neurons. 

On what concerns to the phenomenon of SR itself, an important aspect that should be emphasized is the fact 
that the model we have presented exhibits SR in both monostable and bistable situations, depending on the values 
of the parameters. This remarkable feature contrasts with previous studies for which SR has been found only for 
monostable or bistable systems. Moreover, some bistable systems undergoing SR may become monostable, but under 
this circumstance SR does not take place. In fact, the model under consideration describes a system with a threshold, 
accounting for the firing of the neurons that may behave as a monostable or a bistable system. Therefore, we 
have envisaged a model that may exhibit three different situations (monostable, bistable, and threshold) where the 
phenomenon occurs. 
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FIG. 1. (a) Representation of the potential function U(V) [Eq. ((?])] for = 2 (continuous line), = 3 (dotted line), 
and = 1 (dashed line). The values of the remaining parameters are C = 1, v = 10, 9 = 2, e = 2. (b) Representation of 
U{V) for 9 = 2 (continuous line), 9 = 1.5 (dotted line), and 9 = 2.5 (dashed line). The values of the remaining parameters are 
C = l, u = W, R' 1 = 2, e = 2. 
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FIG. 2. (a) SNR corresponding to Eq. (^) as a function of the noise level for tuo/2n = 0.1 (circles) and luq/2tt = 0.01 
(triangles). The values of the remaining parameters are k = 2, a = 0.5, e = 2, xo = 2, and v — 10. (b) SNR corresponding to 
Eq. (Ua) as a function of the noise level for ujq/2-k = 0.1 (circles) and ujo/2-k = 0.01 (triangles). The values of the remaining 
parameters are k = 2, a — 0.5, e = 2, xq = 2, and f = 10. In both cases the solid line is a guide for the eye. 
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FIG. 3. (a) Representation of the potential function U(V) [Eq. (|jj)]. Same situation as in Fig. |l|a, but e = 8. (b) 
Representation of f7(V). Same situation as in Fig. |l|b, but e = 8. 
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FIG. 5. (a) Time series corresponding to Eq. @ for D = 0.56 (top), D = 5.6 (middle), D = 100 (bottom). The values 
of the remaining parameters are u)q/2h = 0.01, n — 2, a = 0.5, e = 8, xq — 2, and ^ = 10. (b) Time series corresponding to 
Eq. ( |45| ) for D = 0.56 (top), D = 3.3 (middle), D = 100 (bottom). The values of the remaining parameters are ujq/2ty = 0.01, 
K = 2, a — 0.5, e = 8, xq = 2, and v — 10. 
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